Chapter 4 

Functional Field Integral 



In this chapter, the concept of path integration is generalized to integration over quan- 
tum fields. Specifically we will develop an approach to quantum field theory that takes as 

its starting point an integration over all confi.gurations of a given field, weighted by an 
appropriate action. To emphasize the importance of the formulation which, methodologic- 
cally, represents the backbone of the remainder of the text, we have pruned the discussion 
to focus only on the essential elements. This being so, conceptual aspects stand in the 
foreground and the discussion of applications is postponed to the following chapters. 

In this chapter, the concept of path integration will be extended from quantum me- 
chanics to quantum field theory. Our starting point will be from a situation very much 
analogous to that outlined at the beginning of the previous chapter. Just as there are 
two different approaches to quantum mechanics, quantum field theory can also be formu- 
lated in two different ways; the formalism of canonically quantised field operators, and 
functional integration. As for the former, although much of the technology needed to effi- 
ciently implement this framework — essentially Feynman diagrams — originated in high 
energy physics, it was with the development of condensed matter physics through the 50s, 
60s and 70s that this approach was driven to unprecedented sophistication. The reason 
is that, almost as a rule, problems in condensed matter investigated at that time neces- 
sitated perturbative summations to infinite order in the non-trivial content of the theory 
(typically interactions) . This requirement led to the development of advanced techniques 
to sum (subsets of) the perturbation series in many-body interaction operators to infinite 
order. 

In the 70s, however, essentially non-perturbative problems began to attract more and 
more attention — a still prevailing trend — and it turned out that the formalism of 
canonically quantised operators was not tailored to this type of physics. By contrast, the 
alternative approach to many-body problems, functional integration, is ideally suited! 
The situation is similar to the one described in the last chapter where we saw that 
the Feynman path integral provided an entire spectrum of novel routes to approach- 
ing quantum mechanical problems (controlled semi-classical limits, analogies to classical 
mechanics, statistical mechanics, concepts of topology and geometry, etc.). Similarly, the 
introduction of functional field integration into many-body physics spawned plenty of 
new theoretical developments, many of which were manifestly non-perturbative. More- 
over, the advantages of the path integral approach in many-body physics is even more 
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pronounced than in single particle quantum mechanics. Higher dimensionality introduces 
fields of a more complex internal structure allowing for non-trivial topology while, at the 
same time, the connections to classical statistical mechanics play a much more important 
role than in single particle quantum mechanics. 
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Figure 4.1: The concept of field integration. Upper panels: path integral of quantum 
mechanics — integration over all time-dependent configurations of a point particle degree 
of freedom leads to integrals over curves. Lower panels: field integral — integration 
over time dependent configurations of (i-dimensional continuum mappings (fields) leads 
to integrals over generalized {d + 1) -dimensional surfaces. 

All of these concepts will begin to play a role in subsequent chapters when applications 
of the field integral are discussed. Before embarking on the quantitative construction — 
the subject of the following sections — let us first anticipate the kind of structures that one 
should expect. In quantum mechanics, we were starting from a single point particle degree 
of freedom, characterized by some coordinate q (or some other quantum numbers for that 
matter). Path integration then meant integration over all time-dependent configurations 
q(t), i.e. a set of curves t h-h> q(t) (see Fig. 4.1 upper panel). By contrast, the degrees of 
freedom of field theory are continuous objects $(2;) by themselves, where x parameterizes 
some (i-dimensional base manifold and $ takes values in some target manifold (Fig. 4.1, 
lower panel). The natural generalization of a 'path' integral then implies integration over 
a single copy of these objects at each instant of time, i.e. we shall have to integrate 
over generalized surfaces, mappings from {d + l)-dimensional space-time into the field 
manifold, (x, t) 1-^ ^(x, t). While this notion may sound worrying, it is important to 
realize that, conceptually, nothing much changes in comparison with the path integral: 
instead of a one-dimensional manifold — a curve — our object of integration will be a 
{d + l)-dimensional manifold. 

We now proceed to formulate these ideas in quantitative terms. 

> Exercise. If necessary, recapitulate the general construction scheme of path integrals 
(section 3.2.3) and the connection between quantum fields and second quantized operators. 
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4.1 Construction of the Many— body Path Integral 

The construction of a path integral for field operators follows the general scheme outlined 
at the end of section 3.2.3. The basic idea is to segment the time evolution of a quan- 
tum (many-body) Hamiltonian into infinitesimal time slices and to absorb as much as is 
possible of the quantum dynamical phase accumulated during the short time propagation 
into a set of suitably chosen eigenstatcs. But how should these eigenstatcs be chosen? In 
the context of single particle quantum mechanics, the basic structure of the Hamiltonian 
suggested the choice of a representation in terms of coordinate and momentum eigen- 
states. Now, given that many particle Hamiltonians are conveniently expressed in terms 
of creation/ annihilation operators, an obvious idea would be to search for eigenstates of 
these operators. Such states indeed exist and are called coherent states. 



4.1.1 Coherent States (Bosons) 



Our goal is, therefore, to find eigenstates of the Fock space (non-Hermitian) operators 
and a. Although the general form of these states will turn out to be the same for 
bosons and fermions, there are major differences regarding their algebraic structure. The 
point is that the anticommutation relations of fermions require that the eigenvalues of 
an annihilation operator themselves anticommute, i.e. they cannot be ordinary numbers. 
Postponing the introduction of the unfamiliar concept of anticommuting 'numbers' to the 
next section, we first concentrate on the bosonic case where problems of this kind do not 
arise. 

So what form do the eigenstates \4>) of the bosonic Fock space operators a, and 
take? Being a state of the Fock space, an eigenstate |0) can be expanded as 



10)= 



ni,n2, 



|ni,n2---), |ni,n2---) = ^^^^---lO), 



ni,n2, 



where aj creates a boson in state i, Cni,n2,--- represents a set of expansion coefficients, 
and |0) represents the vacuum. Here, for reasons of clarity, it is convenient to adopt this 
convention for the vacuum as opposed to the notation 1^2) used previously. Furthermore, 
the many-body state |ni,n2 • • •) is indexed by a set of occupation numbers: rii in state 
|1), n2 in state |2), and so on. Importantly, the state |0) can, in principle (and will in 
practice) contain a superposition of basis states which have different numbers of particles. 
Now, if the minimum number of particles in state |0) is no, the minimum of a|^|0) must 
be no + 1: Clearly the creation operators a| themselves cannot possess eigenstates. 

However, with annihilation operators this problem does not arise. Indeed, the annihi- 
lation operators do possess eigenstates known as bosonic coherent states 




(4.1) 



where the elements of = {4>i} represent a set of complex numbers. The states \4>) are 
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eigenstates in the sense that, for all i, 



(4.2) 



i.e. they simultaneously diagonalise all annihilation operators. Noting that and a], with 
j 7^ i, commute, Eq. (4.2) can be verified by showing that aexp(0a 



iexp(0a 



Although not crucial to the practice of functional field integration, in the construction of 
the many-body path integral, it will be useful to assimilate some further properties of 
coherent states. 

> By taking the Hermitian conjugate of Eq. (4.2), we find that the 'bra' associated 
with the 'kef |0) is a left eigenstate of the set of creation operators, i.e. for all i, 

(014 = 



(4.3) 



where (f)i is the complex conjugate of (pi, and 



expEi Mi]- 



t> It is a straightforward matter — e.g. by a Taylor expansion of Eq. (4.1) — to show 
that the action of a creation operator on a coherent state yields the identity 



(4.4) 



Reassuringly, it may be confirmed that Eqs. (4.4) and (4.2) are consistent with the 
commutation relations [0^,0]] = Sij: [aj,a]]|0) = (90. 0^ - (j)id^.)\(f)) = Sij\(j)). 

> Making use of the relation (6'|0) = (O|e^'^^"^|0) = e^i^*''^'(O|0) one finds that the 
overlap between two coherent states is given by 



(^10) = exp 



(4.5) 



> Prom this result, one can infer that the norm of a coherent state is given by 



(010) = exp 



y^0i0» 



(4.6) 



> Most importantly, the coherent states form a complete — in fact an overcomplete 
— set of states in Fock space: 



(4.7) 



where (i0jd0i = dRe0i(flm0j, and Ijf represents the unit operator or identity in the 

Fock s]:)ace. 




^Using the result [a, (a'^)"] = n(a^)" ^ (cf. Eq. ??) a Taylor expansion shows aexp{(j)a^)\0) = 



[a, exp (0a 
0exp((/)a1')|O). 



Er=i^(«^r-'io) 
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> Info. The proof of Eq. (4.7) proceeds by straightforward appHcation of Schur's lemma 
(cf. our discussion of the completeness of the spin coherent states in the previous chapter): 
The operator family {a^}, {a|} acts irreducibly in Fock space. According to Schur's lemma, the 
proportionality of the left hand side of Eq. (4.7) to the unit operator is, therefore, equivalent to 
its commutativity with all creation and annihilation operators. Indeed, this property is easily 
confirmed: 

a, J d{l<P)e-^^^^'t>^\<P){<P\ = J d(^,0)e-^^^*'^^</>,|0)(0| = -J d{4>A) (%e-^*^*^') \4>){ct>\ 
2^^^ j (5^^(<^|) = j d{4>,ct>)e-^^^^'^^<t>){<t>\ai, (4.8) 

where, for brevity, we have set (i(0, (j)) = Y\ - dcpidcpi/Tr. Taking the adjoint of Eq. (4.8), one may 
further check that the left hand side of (4.7) commutes with the set of creation operators, i.e. 
it must be proportional to the unit operator. To fix the constant of proportionality, one may 
simply take the overlap with the vacuum: 

J d{4>,(p)e-^i^*'^*{0\4>){4>\0) = J d{4>,(p)e-^i^''^' = 1, (4.9) 

where the last equality follows from Eq. (3.11). Taken together, Eqs. (4.8) and (4.9) prove 
(4.7). Note that the coherent states are overcomplete in the sense that they are not mutually 
orthogonal (sec Eq. (4.5)). The exponential weight e^^^'^''^' appearing in the resolution of the 
identity compensates for the overcounting achieved by integrating over the whole set of coherent 
states. 



With these definitions we have all that we need to construct the many-body path inte- 
gral for the bosonic system. However, before doing so, we will first introduce the fermionic 
version of the coherent state. This will allow us to construct the path integrals for bosons 
and fermions simultaneously, thereby emphasising the similarity of their structure. 



4.1.2 Coherent States (Fermions) 



Surprisingly, much of the formalism above generalises to the fermionic case: As before, it is 
evident that creation operators cannot possess eigenstates. Following the bosonic system, 
let us suppose that the annihilation operators are characterised by a set of coherent states 
such that, for all i, 



ai\v) = Vi\v) 



(4.10) 



where rji is the eigenvalue. Although the structure of this equation appears to be equivalent 
to its bosonic counterpart (4.2) it has one frustrating feature: Anticommutativity of the 
fermionic operators, [a,, aj]+ = 0, where i ^ j, implies that the eigenvalues rji also have 
to anticommute, 



ViVj 



-VjVi 



(4.ii: 



Clearly, these objects cannot be ordinary numbers. In order to define a fermionic version 
of coherent states, we now have two choices: We may (a) accept Eq.(4.11) as a working 
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definition and pragmatically explore its consequences, or (b), first try to remove any 
mystery from the definitions (4.10) and (4.11). This latter task is tackled in the info 
block below where objects {rji} with the desired properties are defined in a mathematically 
consistent manner. Readers wishing to proceed in a maximally streamlined manner may 
skip this exposition and directly turn to the more praxis-oriented discussion below. 

> Info. There is a mathematical structure ideally suited to generalize the concept of 
ordinary number(fields), namely algebras. An algebra .4 is a vector space endowed with a 
multiplication rule Ax A. So, let us construct an algebra A by starting out from a set of 
elements, or generators, r)i e A,i = 1, . . . N, and imposing the rules: 

(i) The elements rji can be added and multiplied by complex numbers, viz. 

Co + CiTji + CjTjj , e A co,Ci,CjeC, (4.12) 
i.e. v4 is a complex vectorspace. 

(ii) The product, A x A ~> A,{rii,rij) ^-^ ViVji associative and anticommutative, i.e. it 
obeys the anti-commutation relation (4.11). Because of the associativity of this operation, 
there is no ambiguity when it comes to forming products of higher order, i.e. (??ir/j)?7fc = 
iJiiVj'nk) = ViVjVk- The definition requires that products of odd order in the number 
of generators anti-commute, while (even,even) and (even,odd) combinations commute 
(exercise) . 

By virtue of (i) and (ii), the set A of all linear combinations cq+X^^i Ylii ■■■in=i ^i,-,inVii ■ ■ - Vin^ 

CQ.Ci-^ G C spans a finite-dimensional associative algebra A,'^ known as the Grsissmann 

algebra'^ (and sometimes also the exterior algebra). 

For completeness we mention that Grassmann algebras find a number of realizations in math- 
ematics, the most basic being exterior multiplication in linear algabra: Given an A/"-dimensional 
vector space V, let V* be the dual space, i.e. the space of all linear mappings, or 'forms' 
A : V ^ C,v A(v), where v ^ V. (Like V, V* is a vector space of dimension N.) Next, define 
exterior multiplication through, (A, A') ^ A A A', where A A A' is the mapping 

AAA' -.V xV C 

{v,v') ^ A{v)A'{v')-A{v')A'{v). 

This operation is manifestly anti-commutative: A A A' = —AAA'. Identifying the N linear basis 
forms Aj ^ iji with generators and A with the product, we see that the space of exterior forms 
of a vector space forms a Grassmann algebra. 



whose dimension can be shown to be 2^ (exercise) 



Hermann Giinter Grassmann 

1809-1877: credited for inventing 
what is now called Exterior Alge- 
bra. 
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Apart from their anomalous commutation properties, the generators {r)i}, and their 
product generahzations {f]if]j, r]ir]jr]k, ■ ■ .} resemble ordinary, albeit ant i- commutative num- 
bers. (In practice, the algebraic structure underlying their definition can safely be ignored. 
All we will need to work with these objects is the basic rule (4.11) and the property (4.12).) 
We emphasize that A not only contains anticommuting but also commuting elements, i.e. 
linear combinations of an even number of Grassmann numbers r)i are overall commutative. 
(This mimics the behaviour of the Fock space algebra: products of an even number of 
annihilation operators ajflj . . . commute with all other linear combinations of operators 
Gi. In spite of this similarity, the 'numbers' rji must not be confused with the Fock space 
operators; there is nothing on which they act.) 

To make practical use of the new concept, we need to go beyond the level of pure 
arithmetic. Specifically, we need to introduce functions of anti-commuting numbers, and 
elements of calculus. Remarkably, most of the concepts of calculus not only naturally 
generalize to anti-commuting number fields, but contrary to what one might expect, the 
anti-commutative generalization of differentiation, integration, etc. turns out to be much 
simpler than in ordinary calculus. 



> Functions of Grassmann numbers are defined via their Taylor expansion: 

k 

5=0 



/(6,---,e.) = E E -^ Qc Qc. ,_a.-..en, (4.13) 
n=0 n,-in=l ■ "^^i ■ ■ ■ 



where / is an analytic function. Note that the anticommutation properties of the 
algebra implies that the series terminates after a finite number of terms. For exam- 
ple, in the simple case where rj is first order in the generators of the algebra, N — 1, 
and f{ri) = /(O) + f{Q)ri (since 77^ = 0). 

> Differentiation with respect to Grassmann numbers is defined by 

dmVj = Sij (4.14) 



Note that in order to be consistent with the commutation relations, the differential 
operator 5^. must itself be anti-commutative. In particular, druTjjTji —rjj. 

> Integration over Grassmann variables is defined by 



J drii = 0, J driiTii = 1 



(4.15) 



Note that the definitions (4.13), (4.14) and (4.15) imply that the action of Grass- 
mann differentiation and integration are effectively identical, viz. 

/ dvfiv) = J dvifiO) + nO)rj) = /'(O) = dj{rj) . 
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With this background, let us now proceed to apply the Grassmann algebra to the 
construction of fermion coherent states. To this end we need to enlarge the algebra even 
further so as to allow for a multiplication of Grassmann numbers by fermion operators. 
In order to be consistent with the anticommutation relations, we need to require that 
fermion operators and Grassmann generators anticommute. 



0. 



(4.16) 



It then becomes a straightforward matter to demonstrate that fermionic coherent 
states are defined by 



\v) = exp 




|0) 




i 





(4.17) 



i.e. by a structure perfectly analogous to the bosonic states (4.1)."^ It is a straightforward 
matter — and also a good exercise — to demonstrate that the properties (4.3), (4.4), 
(4.5), (4.6) and, most importantly, (4.7) carry over to the fermionic case. One merely has 
to identify aj with a fermionic operator and replace the complex variables 0j by rji G A. 
Apart from a few sign changes and the ^-valued arguments, the fermionic coherent states 
differ only in two respects from their bosonic counterpart: firstly, the Grassmann variables 
fji appearing in the adjoint of a fermion coherent state. 



(^1 



exp 



(0| exp 



are not related to the rjiS of the state \ri) via some kind of complex conjugation. Rather rji 
and fji are strictly independent variables.^ Secondly, the Grassmann version of a Gaussian 
integral (exercise), J dfjdr]e~^ — 1 does not contain the factors of tt characteristic of 
standard Gaussian integrals. Thus, the measure of the fermionic analogue of Eq. (4.7) 
does not contain a tt in the denominator. 

For the sake of future reference, the most important properties of Fock space coherent 
states are summarised in table 4.1. 



""^To prove that the states (4.17) indeed fulfil the defining relation (4.10), we note that 



aiexp(-77ia|)|0) ^ = ^ ai{l - r,ial)\0) ^ = ' V^a^al\0) = V^\0) = V^i'^ " V^alM = r?i exp(-r?ia^)|0). This, 
in combination with the fact that Oj and r/jfij {i ^ j) commute proves (4.10). Note that the proof has 
actually been simpler than in the bosonic case. The fermionic Taylor series terminates after the first 
contribution. This observation is representative of a general rule: Grassmann calculus is simpler than 
standard calculus; all series are finite, integrals always converge, etc. 

'''In the literature, complex c;onjugation of Grassmann variables is sometimes defined. Although ap- 
pealing from an aesthetic point of view — symmetry between bosons and fermions — this concept is 
problematic. The difficulties become apparent when supersymmetric theories are considered, i.e. 
theories where operator algebras contain both bosons and fermions (tlic^ so c;alled super algebras). It is 
not possible to introduce a complex conjugation that leads to compatibility with the commutation rela- 
tions of a super-algebra. It therefore seems to be better to abandon the concept of Grassmann complex 
conjugation altogether. Note that although, in the bosonic case, complex conjugation is inevitable (in 
order to define convergent Gaussian integrals, say), no such need arises in the fermionic case. 
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Definition 


|V') = exp Cj^M^ |0) 

i 


Action of Qi 


ai\ijj) ^i/jilijj), {'il}\ai ^ d^.{il)\ 


Action of a| 


al|^) = Ca^jV'), {^^J\a\ = {^P\i,, 


Overlap 


= exp ^^l^i 

Li J 


Completeness 


1 d{ij,ij) e-^i^'^^\ij){^\ = U 



Table 4.1: Basic properties of coherent states for bosons {( = 1, ipi E C) and fermions 
(C = —1, ipi e A). In the last line, the integration measure is defined as dl'ipjip) = 

i ^(i+0/2 ■ 

> Info. Grsissmann Gaussian Integration: Finally, before turning to the development 
of the functional field integral, it is useful to digress on the generalization of higher dimen- 
sional Gaussian integrals for Grassmann variables. The prototype of all Grassmann Gaussian 
integration formulae is the identity 

(4.18) 

where a G C takes arbitrary values. Eq. (4.18) is derived by a first order Taylor expansion of 
the exponential and application of Eq. (4.15). The multi-dimensional generalization of (4.18) is 
given by 

' d{^, 0)e-^^^'^ = det A , (4.19) 




where (f) and (f) are A^-component vectors of Grassmann variables, the measure d{(f), (p) = 
YliLi d4>id(j)i, and A may be an arbitrary complex matrix. For matrices that are unitarily diag- 
onalisable, A = U^DU, with U unitary, and D diagonal, Eq. (4.19) is proven in the same way 
as its complex counterpart (3.17): One changes variables (f) V^(t>, (p ~^ U"^(^. Since detU = 1, 
the transform leaves the measure invariant (see below) and leaves us with N decoupled integrals 
of the type (4.18). The resulting product of N eigenvalues is just the determinant of A (cf. 
the later discussion of the partition function of the non-interacting gas). For general (non- 
diagonalisable) A, the identity is established by a straightforward expansion of the exponent. 
The expansion terminates at A'^th order and, by commuting through integration variables, it may 
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be shown that the resulting A'^th order polynomial of matrix elements of A is the determinant.^ 
Keeping the analogy with ordinary commuting variables, the Grassmann version of Eq. (3.18) 

reads 



e^^^ ''^detA 



(4.22) 



To prove the latter, we note that J drjf^r]) = J dij f (rj + u) , i.e. in Grassmann integration one can 
shift variables as in the ordinary case. The proof of the Gaussian relation above thus proceeds 
in complete analogy to the complex case. As with Eq. (3.18), Eq. (4.22) can also be employed 
to generate further integration formulae. Defining (• • •) = detA"^/ d{^,(f))e'~'^ '^'^(" " ")> 
expanding both the left and the right hand side of (4.22) to leading order in the 'monomial' 
UjUi, one obtains {rijfji) = Aj^. Finally, the A'^-fold iteration of this procedure gives 



{VjiVj2 ■ ■ ■ VjnVhVi2 ■ ■ ■ = sgn P) Aj^]^^ . . . A.\^^ 



where the signum of the permutation accounts for the sign changes accompanying the interchange 
of Grassmann variables. Finally, as with Gaussian integration over commuting variables, by 
taking N ^ oo, the Grassmann integration can be translated to a Gaussian functional integral. 



4.2 Field Integral for the Quantum Partition Func- 
tion 

Having introduced the coherent states, we will see that the construction of path inte- 
grals for many-body systems no longer presents substantial difficulties. However, be- 

^As with ordinary integrals, Grassmann integrals can also be subjected to variable transforms. 
Suppose we are given an integral / d{<l>, 4>)f{4>i 4') and wish to change variables according to 

z/ = M^, z/ = M>, (4.20) 

where, for siniplic;ity, M and M' are complex matrices (i.e. we here restrict ourselves to linear transforms). 
One can show that 

Pi . . . z/jv = (det M)(^i . . . ^jv, i/i . . . z/jv = (det M')0i . . . 0jv • (4-21) 

(There are different ways to prove this identity. The most straightforward is by explicitly expanding (4.20) 
in components and commuting all Grassmann variables to the right. A more elegant way is to argue that 
the coefficient relating the right and the left hand sides of (4.21) must be an A'^th order polynomial of 
matrix elements of M. In order to be consistent with the aiiti c;onnmitation behaviour of Grassmann 
variables, the polynomial must obey commutation relations which uniquely characterise a determinant. 
Excercise: Check the relation for N = 2.) On the other hand, the integral of the new variables must obey 
the defining relation, J dvvi . . . Pn = J dvvi . . . vm = {~)^'^^ , where dP = Ilj=i '^'^i the sign on the 
right hand side is attributed to ordering of the integrand, viz. / di>\dv2ViV2 = — J dvivi j dv2V2 = — 1- 
Together Eqs. (4.21) and (4.20) enforce the identities dv = {dei M)' ^ dcj), dv = (det M')"^#, which 
combine to give 

j d(,^,(/.)/(,^,</.)=det(MM') j d{v,i,')f{4>{v),4>{v)). 
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fore proceeding, we should address the question; what does the phrase 'path-integral for 
many-body systems' actually mean? In the next chapter we will see that much of the 
information about a quantum many-particle systems is encoded in expectation values of 
products of creation and annihilation operators, i.e. expressions of the structure {a^a . . .). 
By an analogy to be explained then, objects of this type are generally called correlation 
functions. More important for our present discussion, at any finite temperature, the 
average (. . .) entering the definition of the correlation function runs over the quantum 
Gibbs^ distribution p = e~^''^~^^'> / Z, where, as usual, 

Z = tr e-^^^-'^^) = Y.{n\e-^^^-^^^ \n) , (4.23) 

n 

is the quantum partition function, j3 = 1/T, /i denotes the chemical potential, and the 
sum extends over a complete set of Fock space states (For the time being we neither 

specify the statistics of the system — bosonic or fermionic — nor the structure of the 
Hamiltonian.) 

Ultimately, we will want to construct and evaluate the path integral representations of 
many-body correlation functions. Later wc will see that all of these representations can be 
derived by a few straightforward manipulations form a prototypical path integral, namely 
that for Z. Further, the (path integral of the) partition function is of importance in its 
own right: It contains much of the information needed to characterise the thermodynamic 
properties of a many-body quantum system.^ We thus begin our journey into many-body 
field theory with a construction of the path integral for Z. 

To prepare the representation of the partition function (4.23) in terms of coherent 
states, one must insert the resolution of identity 

Z^ I d{i),i))e-^i^''^'^{n\i)){'^\e-^^"-''^^\n) . (4.24) 

We now wish to get rid of the — now redundant — Fock space summation over \n) 
(another resolution of identity). To bring the summation to the form = 1:f, 



Josiah Willard Gibbs 1839 1903: cred- 
ited with the development of chemical 
thermodynamics, he introduced concepts 

of free energy and chemical potential. 




^In fact, the statement above is not entirely correct. Strictly speaking thermodynamic properties 
involve the thermodynamic potential Q. = —T In Z rather than the partition function itself. At first 
sight it seems that the difference between the two is artificial one might first calculate Z and then 
take the logarithm. However, typically, one is unable to determine Z in closed form, but rather one 
has to perform a perturbative expansion, i.e. the result of a calculation of Z will take the form of a 
series in some small parameter c. Now a problem arises when the logarithm of the series is taken. In 
particular, the Taylor series expansion of 2^ to a given order in e does not automatically determine the 
expansion of fl to the same order. Fortunately, the situation is not all that bad. It turns out that the 
logarithm essentially rearranges the combinatorial structure of the perturbation series in an order known 
as a cumulant expansion. 
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one must commute the factor to the right hand side. However, in performing 

this seemingly innocuous operation, we must be careful not to miss a sign change whose 
presence will have important consequences for the structure of the fermionic path integral: 
Indeed, it may be checked that, whilst for bosons, {n\%l)) {■il)\n) = (ipln) (nlip) , the fermionic 
coherent states change sign upon permutation, {n\ip){ip\n) = {—tp\n) {n\tp) (i.e. {—ip\ = 
exp (— i^iO'i))- The presence of the sign is a direct consequence of the anti-commutation 
relations between Grassmann variables and Fock space operators (exercise). Note that, 
as both H and contain elements even in the creation/annihilation operators, the sign 
is insensitive to the presence of the Boltzmann factor in (4.24). Making use of the sign 
factor C, the result of the interchange can be formulated as the general expression 

Z = /"d(V;,V')e-^*'^^^*^(CV'|e-^(^-'^'^)|n)(n|V') 

n 

d{i^,i^) e-^i'^'''^'{C,i^\e-^^"-''^^\i^) , (4.25) 



where the equality is based on the identity \n) {n\ = Ijc. Eq. (4.25) can now be directly 
subjected to the general construction scheme of the path integral. 

To be concrete, let us assume that the Hamiltonian is limited to a maximum of two- 
body interactions (cf. Eqs. (2.5) and (2.9)), 

H{a\ a) = ^ hija\a^ + ^ Vijki ■ (4-26) 

ij ijkl 

Note that, to facilitate the construction of the field integral, it is helpful to arrange for all 
of the annihilation operators to stand to the right of the creation operators. Fock space 
operators of this structure are said to be normal ordered.^ The reason for emphasising 
normal ordering is that such an operator can be readily diagonalised by means of coher- 
ent states: Dividing the 'time interval' j3 into segments and inserting coherent state 
resolutions of identity (steps 1, 2 and 3 of the general scheme), Eq. (4.25) assumes the 
form 



r ^ 

^ j\r n 



y,o=f^JV n=0 



simi- 



where 5 = f3/N and ^^^^^g^ = E., /^..^.^j + V^ki A^.^ = H{^, 
larly N^ipjip')). Here, in writing Eq. (4.27), we have adopted the shorthand ip^ = {'?/'"}, 
etc. Finally, sending N ^ oo and taking limits analogous to those leading from (3.5) to 



^More generally, an operator is defined to be 'normal ordered' with respect to a given vacuum state 
|0), if and only if, it annihilates |0). Note that the vacuum need not necessarily be defined as a zero 
particle state. If the vacuum contains partic;l(^s. normal ordering need not lead to a representation where 
all annihilators stand to the right. If, for whatever reason, one is given a Hamiltonian whose structure 
differs from (4.26), one can always affect a normal ordered form at the expense of introducing commutator 
terms. For example, normal ordering the quartic term leads to the appearance of a quadratic contribution 
which can be absorbed into hai3- 
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(3.6) we obtain the continuum version of the path integral, 



10 



(4.28) 



where D{il)^ ip) = hniAr^oo 11^=1 d{i/j^, ip"'), and the fields satisfy the boundary condition 

m^Ci^iP), m-W)- (4-29) 

Written in a more explicit form, the action associated with the general pair-interaction 
Hamiltonian (4.26) can be cast in the form 



dr 



. ij ijkl 



(4.30) 



Notice that the structure of the action fits nicely into the general scheme discused 
in the previous chapter. By analogy, one would expect that the exponent of the many- 
body path integral carries the significance of the Hamiltonian action, S ~ j {p<i — H), 
where {q,p) symbolically stands for a set of generalized coordinates and momenta. In 
the present case, the natural pair of canonically conjugate operators is (a, a^). One would 
then interpret the eigenvalues (■0, ■0) as 'coordinates' (much as {q,p) are the eigenvalues of 
the operators {q,p)). Adopting this interpretation, we see that the exponent of the path 
integral indeed has the canonical form of a Hamiltonian action and, therefore, is easy to 



memorize. 



Eqs. (4.28) and (4.30) define the functional integral in the time representation 

(in the sense that the fields are functions of a time variable) . In practice we shall mostly 
find it useful to represent the action in an alternative, Fourier conjugate representation. 
To this end, note that, due to the boundary conditions (4.29), the functions iP{t) can be 
interpreted as functions on the entire Euclidean time axis which are periodic/antiperiodic 
on the interval [0, /?]. As such they can be represented in terms of a Fourier series. 



1 



where 



2nnT, 

(2n + l)7rT, 



bosons, 
fermions 



neZ 



(4.31) 



are known as Matsubara frequencies. Substituting this representation into (4.28) 
and (4.30), we obtain 2 = J D{ip,ip)e~^^'^''^\ where D{ip,ip) = Ylndi'^m'^n) defines the 



^"^Whcroas the bosonic contiimuin limit is indeed perfectly equivalent to the one taken in constructing 
the quantum mechanical path integral (lim^^o <5~^(V'"^^ ~ V'") = drlr^nS gives an ordinary derivative 
etc.), a novelty arises in the fermionic case. The notion of replacing differences by derivatives is purely 
symbolic for Grassmann variables. There is no sense in which -0"+^ — t/j" is small. The symbol dr'tp rather 
denotes the formal (and well defined expression) lim^^o 6~^{'tp"'^^ — xp^). 
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measure (for each Matsubara index n we have an integration over a coherent state basis 
{iV'n)})/^ and the action takes the form 

"0] = XI [{-'^^n - lA Sij + hij] Ipjn + X Vijkl 1pini1pjn2'ipkn3'iplnJni+n2,n3+ni4:.S2) 

Here we have used the identity JJ^ (ire""^"'^ = P^u>n,o- Eq. (4.32) defines the frequency 
representation of the action. 

> Info. In performing calculations in the Matsubara representation, one sometimes runs 
into convergence problems (which will manifest themselves in the form of ill-convergent Matsub- 
ara frequency summations): In such cases it will be important to remember that Eq. (4.32) does 
not actually represent the precise form of the action. What is missing is a convergence generat- 
ing factor whose presence follows from the way in which the integral was constructed, and which 
will save us in cases of non-convergent sums (except, of course, in cases where divergences have 
a physical origin). More precisely, since the fields -ijj are evaluated infinitesimally later than the 
operators ■0 (cf. Eq. (4.27)), the h and /i dependent contributions to the action acquire a factor 
exp(— icJnJ), 6 infinitesimal. Similarly, the V contribution acquires a factor exp{—i{ujni +^n2)^)- 
In cases where the convergence is not critical, we will omit these contributions. However, once 
in a while it is necessary to remember their presence. 



4.2.1 Partition Function of Non— Interacting Gas 

As a first exercise, let us consider the quantum partition function of the non-interacting 
gas. (Later, this object will prove useful as a 'reference' in the development of weakly 
interacting theories.) In some sense, the field integral formulation of the non-interacting 
partition function has a status similar to that of the path integral for the quantum har- 
monic osciallator: The direct quantum mechanical solution of the problem is straightfor- 
ward and application of the full artillery of the field integral seems somewhat ludicrous. 
Prom a pedagogical point of view, however, the free partition function is a good problem; 
it provides us with the welcome opportunity to introduce a number of practical concepts 
of field integration within a comparatively simple setting. Moreover, the field integral 
representation of the free partition function will be an important operational building 
block for our subsequent analysis of interacting problems. 

Consider, then the partition function (4.28) with Ho{tjj,ijj) = Yliij^iHo,ij'4'j- Diag- 
onalising Hq by a unitary transformation [/, Hq — UDW and transforming integration 
variables Wip = (f), the action assumes the form, S = 4>an{—'ii^n + Co) (pan, where 

= f„ — /' and f„ arc the single particle eigenvalues. Remembering that the fields (pai'T) 

^^Notice, however, that the fields ^/^n carry dimension [energy]"^/^, i.e. the frequency coherent state 
integral is normalized as / d{'ipn,'4'n) e"'^"'^'^" = (/3e)~^. 

^^As for the signs of the Matsubara indices appearing in Eq. (4.32), note that the Fourier representation 
of i/j is defined as V'(r) = V'^e"''"^"'^, V'n = dr 'ip{T)e~^'^'^ . In the bosonic case, this sign 

convention is motivated by -(/i being the complex conjugate of V'- For reasons of notational symmetry, this 
convention is also adopted in the fermionic case. 
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are independent integration variables (exercise: why does the transformation ^0 — > have 
a Jacobian of unity?), we find that the partition function decouples, Z — Yla-^a, where 

Z„ = /" D{^,, e-^n^'^-(-'-n+ia)4>an ^ ^[[^{-iuj^ + ^Jj'C , (4.33) 

n 

and the last equality follows from the fact that the integrals over are onc-dimcnsional 
complex or Grassmann Gaussian integrals. Here, let us recall our convention defining 
( = 1(— 1) for bosonic (fermionic) fields. At this stage, we have left all aspects of field 
integration behind us and reduced the problem to one of computing an infinite product 
over factors icUn — Ca- Since products are usually more difficult to get under control than 
sums, we take the logarithm of Z to obtain the free energy 

F = -TlnZ = -rC^ln[/3(-zw„ + Q] . (4.34) 

an 



> Info. Before proceeding with this expression, let us take a second look at the interme- 
diate identity (4.33). Our calcuation showed the partition function to be the product over all 
eigenvalues of the operator —iu + H — defining the action of the non-interacting system 
(here, lo = {con^nn'})- A.S such, it can be written compactly as: 



Z = det 



P{-iu + H- fxN) 



-C 



(4.35) 



This result was derived by first converting to an eigenvalue integration and then performing tlie 
one-dimensional integrals over 'eigencomponents' While technically straightforward, that 
— explicitly representation-dependent — procedure is not well suited to generalization to more 
complex problems. (Keep in mind that later on we will want to embed the free action of the 
non-interacting problem into the more general framework of an interacting theory.) 

Indeed, it is not necessary to refer to an eigenbasis at all: In the bosonic case, Eq. (3.17) 
tells us that Gaussian integration over a bilinear ~ generates the inverse determinant of 
X. Similarly, as we have seen, Gaussian integration extends to the Grassmann case with the 
determinants appearing in the numerator rather than in the denominator (as exemplified by 
(4.35)). (As a matter of fact, (4.33) is already a proof of this relation.) 



We now have to face up to a technical problem: How do we compute Matsubara sums 
of the form X]„ln(ia;„ — x)l Indeed, it takes little imagination to foresee that sums of 
the type Xlni ^i^ni,^n2, ■ ■ whcrc X symbolically stands for some function, will 
be a recurrent structure in the analysis of functional integrals. A good ansatz would 
be to argue that, for sufficiently low temperatures (i.e. temperatures smaller than any 
other characteristic energy scale in the problem), the sum can be traded for an integral, 
viz. TY2^ — > J duj /{2t:). However, this approximation is too crude to capture much of 
the characteristic temperature dependence in which one is usually interested. Yet there 
exists an alternative, and much more accurate way of computing sums over Matsubara 
frequencies: 
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> Info. Consider a single Matsubara frequency summation, 



(4.36) 



where h is some function and ujn may be either bosonic or fermionic (cf. Eq. (4.31)). The 
basic idea behind the standard scheme of evaluating sums of this type is to introduce a complex 
auxiliary function g{z) that has simple poles at z = iiVn- The sum S then emerges as the sum 
of residues obtained by integrating the product gh along a suitably chosen path in the complex 
plane. Typical choices of g include 



exp(/3z) 
/3 



. exp(/5z) + 1 ' 



bosons 



fermions 



and g{z) 



( P 

— coth(/3z/2), bosons 

— tanh.{f3z/2), fermions 
k 2 



(4.37) 



where, in much of this section, we will employ the functions of the first column. (Notice the 
similarity between these functions and the familiar Fermi and Bose distribution functions.) In 
practice, the choice of the counting function is mostly a matter of taste, save for some cases 
where one of the two options is dictated by convergence criteria. 




Figure 4.2: Left: the integration contour employed in calculating the sum (4.36). Right: 
the deformed integration contour. 



Integration over the path shown in the left part of Fig. 4.2 then produces 

^ dzg{z)h{-iz) = C,^ Res (5(2;)/i(-fz))|^^.^^ = ^/i(u;„) = S , 



2TTi 



where, in the third identity, we have used the fact that the 'counting functions' g are chosen so 
as to have residue ( and it is assumed that the integration contour closes at z — > itioo. (The 
difference between using the first and the second column of (4.37) lies in the value of the residue. 
In the latter case, it is equal to unity rather than (^.) Now, the integral along a contour in the 
immediate vicinity of the poles of g is usually intractable. However, as long as we are careful not 
to cross any singularities ol g or the function h{—iz) (the latter symbolically indicated by isolated 
crosses in the figure^^) we are free to distort the integration path, ideally to a contour along 



^■^Remember that a function that is bounded and analytic in the entire complex plane is constant, i.e. 
every 'interesting' function will have singularities. 
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which the integral can be performed. Finding a suitable contour is not always straightforward. 
If the product hg decays sufficiently fast at \z\ oo (i.e. faster than z~^), one will ususally try 
to 'inflate' the original contour to an infinitely large circle (Fig. 4.2, right). ^'^ The integral along 
the outer perimeter of the contour then vanishes and one is left with the integral around the 
singularities of the function h. In the simple case where h{—iz) posesses a number of isolated 
singularities at {zk} (i.e. the situation indicated in the figure) we thus obtain^^ 

S = f H-i^)9{^) = -Cj2 ResM-i%(^)|.=., , (4.38) 

k 

i.e. the task of computing the infinite sum S has been reduced to that of evaluating a finite 
number of residues — a task that is always possible! 

To illustrate the procedure on a simple example, let us consider the function 



where (5 is a positive infinitesimal.^^ To evaluate the sum S = X]„/i(i^n), we first observe that 
the product h{—iz)g{z) has benign convergence properties. Further, the function h{—iz) has a 
simple pole that, in the limit 5 — > 0, lies on the real axis at z = ^. This leads to the result 

^h{ujn) = C'Resg{z)h{-iz)\^^^ = —^pfz^r " 

n ' 

We have thus arrived at the important identity 

(4.39) 




-^"'Notice that the condition linii^i^oo l^ffl < is not as restrictive as it may seem. The reason is 
that the function h will be mostly related to physical observables that approach some limit (or vanish) 
for large excitation energies. This implies vanishing in at least portions of the complex plane. The 
convergence properties of g depend on the concrete choice of the counting function. (Exercise: explore 
the convergence properties of the functions shown in Eq. (4.37).) 

^^If you are confused about signs, note that the contour encircles the singularities of /i in a clockwise 
direction. 

^^Indeed, this choice of h is actually not as artificial as it may seem. The expectation value of the 
number of particles in the grand canonical ensemble is defined through the identity N = —dF/dji 
where F is the free energy. In the non-interacting case, the latter is given by Eq. (4.34) and, remembering 
that = — M) one obtains N w — C^X^q„ . Now, why did we write instead of '='? The 

reason is that the right hand side, obtained by naive differentiation of (4.34), is ill convergent. (The 
sum ^T^' arbitrary, does not exist!) At this point we have to remember the remark made in 

the on page 136, i.e., had we carefully treated the discretisation of the field integral, both the logarithm 
of the free energy and dfj_F would acquire infinitesimal phases cxp(— icij„(5). As an exercise, try to keep 
track of the discretisation of the field integral from its definition to Eq. (4.34) to show that the accurate 
expression for N reads 



an 



where h is the function introduced above. (Note that the necessity to keep track of the lifebuoy e""^"'' 
docs not arise too often. Most Matsubara sums of physical interest relate to functions / that decay faster 
than z~^.) 



Concepts in Theoretical Physics 



140 



CHAPTER 4. FUNCTIONAL FIELD INTEGRAL 



where 



npie) 



1 



"-B(e) 



1 



(4.40) 



g/^Ce-A*) _ 1 



are the Fermi/Bose distribution functions. As a corollary we note that the expectation value 
for the number of particles in a non-interacting quantum gas assumes the familiar form A'^ = 

Before returning to our discussion of the partition function, let us note that life is not always 
as simple as the example above. More often than not, the function h not only contains isolated 
singularities but also cuts or worse singularities. Under such circumstances, finding a good 
choice of the integration contour can be far from straightforward! 



Returning to the problem of computing the sum (4.34), consider for a moment a fixed 
eigenvalue = ct- In this case, we need to evaluate the sum 5* = J2n^i^ri), where 
h{LJn) = CT\ia[P{—iUn + 0] = (T\ia[P{iUn — ^)] + C and C is an inessential constant. As 
discussed before, the sum can be represented as S* = —^<fg{z)h{—iz), where g{z) = 
I3{e^^ _ (^)-i ig (^ times) the distribution function and the contour encircles the poles of 
g as in Fig. 4.2, left. Now, there is an essential difference with the example discussed 
previously, viz. the function h{—iz) = (Tln{z — C,) + C has a branch cut along the real 
axis, z e (— oo,^) (see the figure). To avoid contact with this singularity one must distort 
the integration contour as shown in the figure. Noticing that the (suitably regularized, cf. 
our previous discussion of the particle number N) integral along the perimeter vanishes, 
we conclude that 



where = e±iri, rj is a. positive infinitesimal, and we have used the fact that g{e^) — g{e) 
is continuous across the cut. (Also, without changing the value of the integral (exercise: 
why?), we have enlarged the integration interval from (— cxd,^] to (— cxd, oo)). To evaluate 





)( 
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)( 




Figure 4.3: 
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the integral, we observe that g{e) — —(d^ In (l — (e and integrate by parts: 

^ - -i /* In (1 - Ce-) (j^ - ^) ''^'> CT^ (1 - Ce-«) . 

Insertion of this result into Eq. (4.34) finally obtains the familiar expression 

F = ^ In (1 - (e-^^^'^-'^)) (4.41) 

a 

for the free energy of the non-interacting Fermi/Bose gas. While this result could have 
been obtained much more straightforwardly by the methods of quantum statistical me- 
chanics, we will shortly see how powerful a tool the methods discussed in this section are 
when it comes to the analysis of less elementary problems! 

4.3 Summary and Outlook 

This concludes our preliminary introduction to the field integral. We have learned how to 
represent the partition function of a quantum many-body system in terms of a generalized 
path integral. The field integral representation of the partition function will be the basic 
platform on which all our further developments will be based. In fact, we are now in a 
position to face up to the main problem addressed in this text: practically none of the 
'non-trivial' field integrals in which one might be interested can be executed in closed form. 
This reflects the fact that, save for a few exceptions, interacting many-body problems do 
not admit closed solutions. In the following chapter, we will introduce approximation 
stratgies for addressing interacting theories by exploring some physical applications of 
the field integral. 
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